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Abstract 

X 

Hadron colliders offer a unique opportunity 
to test perturbative QCD because, rather than 
producing events at a specific beam energy, the dynamics 
of the hard scattering is probed simultaneously at a wide 
range of momentum transfers. 
This makes the determination of as 

and the parton density functions (PDF) at hadron colliders particularly interesting. 
In this paper we restrict ourselves to extracting as for a given PDF at a scale 
which is directly related 

to the transverse energy produced in the collision. As an example, we focus on 
the single jet inclusive transverse energy 

distribution and use the published '88-'89 CDF data with an integrated luminosity 
of 4.2 pb -1 . The evolution of the coupling constant over a wide range of scales (from 
30 GeV to 500 GeV) is clearly shown and 



is in agreement with the QCD expectation. The data to be obtained in the current 
Tevatron run (expected to be well in excess 100 pb _1 for both the CDF 
and D0 experiments) will significantly decrease the experimental errors. 
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1 Introduction 



Hadronic collisions at the Fermilab TEVATRON 
offer excellent opportunities to study QCD over a 

broad range of momentum transfers ranging from a few GeV in the transverse 
momentum distribution of the Z boson up to almost half of the beam energy in the single jet 
inclusive transverse energy distribution. While the experiments at LEP and HERA have 
set well defined goals for QCD studies, 

hadron colliders tend to be thought of as discovery machines probing the high energy frontier. For 
example, at Fermilab, the major effort has been concentrated on the study of the Top Quark and VF-mass 
measurements. 

In this paper, we try to redress this imbalance and outline a possible goal for QCD studies at hadron 
colliders. 

Achieving this goal will give both a rigorous 

test of QCD and a reduction of the experimental systematic errors in the other studies at Fermilab. 

One possible goal of QCD studies at the Main Injector 111 is to use the QCD data set to determine the 
input parameters of the theory, 

in other words, as and the parton density functions (PDF's), 

without input from other experiments.^ This should also allow the determination of the gauge symmetry 
responsible for the strong interactions thereby extending similar measurements at LEP Q . 

In order to achieve this goal, we need to make several intermediate steps to identify problems in both 

experiments and theory. The run 1A and IB data can be used to 

gain experience in how to analyze the data and to identify those distributions which can be measured 
accurately and calculated reliably. We will break the program into four steps with each phase contributing 
to a better understanding of QCD at hadron colliders. 

In the first phase we use the PDF's obtained from global analyses ||, ||, |[ |) and the associated as{Mz) 
as input parameters. 

Then by comparing data to theory we 
1 Note that the range of x and Q 2 

probed in hadron-hadron collisions is rather different from that probed at HERA. 
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can identify those cross sections which are most sensitive to the input parameters. 

Using run 1A data it has become clear that certain distributions will be better than others in determining 
the parton density functions and as- For example, the parton density functions can be constrained from 
di-jet data using angular correlations the same-side to opposite-side ratio || ||, and via the triply 
differential cross section |n], 12, [l3| pT, |15| while the strong coupling can be determined from vector boson 
production at large transverse momentum jl^, 17, lSy . 

In the second phase we will assume a given PDF set as being correct and extract as- Measuring as at a 
hadron collider is rather different than measuring as at LEP with the most important difference being the 
fact that one can 

measure as from momentum transfers as low as a few 

GeV all the way up to 500 GeV simultaneously and with high statistics [fl5| . In this paper we will use 

the 

one-jet inclusive transverse energy distribution 

measured from the '88-'89 CDF data with an integrated luminosity of 4.2 pb _1 

to illustrate this method. The analysis can be repeated for the current CDF and D0 data sets 

increasing the integrated luminosity to well over 100 pb _1 . These increased statistics 

will have a major impact on both the statistical and systematic error 

relative to the '88-'89 data set. 

The results in this paper are therefore just 

an illustration of the method and no detailed effort has been made to determine the 

experimental systematic errors thoroughly. This would require detailed knowledge of the 

correlation matrix for the systematic error which is not readily available. While the value of as extracted 



this way cannot be considered on the same footing as that measured at LEP because the PDF's themselves 
are dependent on as-, this measurement will nevertheless provide valuable information. 

For example, the extracted as must be consistent with the as used in the PDF, or else the data is 
incompatible with 

this particular set. If one finds that the extracted as is compatible with the PDF the measurement gives 
an additional constraint on the PDF at large x and Q 2 . Further, one can also study the evolution of as for 
a wide range of momentum transfers. n 



2 An alternative approach has been followed by UA2 [fL6| and is now being pursued actively by both CDF 
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In the third phase we determine both 

the PDF's and as simultaneously using the triply differential di-jet inclusive distributions g, [ll], [T^| , 
possibly including flavor tagging, yielding an a 5 

that is completely independent of the DIS data set. 

In principle the measurement is very simple, 

the parton fractions are determined by summing over the rapidity weighted 
transverse momenta of the particles 
produced by the hard scattering, 

Xi t 2 = (J2i E^e ±rii )/VS |l3|, Q. However, since it is impossible to measure and identify all the 
particles associated with the hard scattering, we are forced to rely on 
higher order calculations to estimate the unobserved 
radiation. It might therefore be prudent to separate 

this phase into two steps by first determining the distribution of gluons in the 

proton and assuming the distribution of charged partons is determined by the DIS data set. The reason 
for this is that in deep inelastic scattering the virtual photon directly probes the charged parton distributions. 
The effects of the gluon distribution enter first 

at next-to-leading order and cause, for example, 

scaling violations in 

the slope of F2. On the other hand, in hadron colliders, 
the gluon density enters at lowest order 

and a more direct measurement should be possible. For example, by using the triply differential di- 
jet data one probes the gluon distribution directly with essentially unlimited statistics. After the gluon 
distribution has been successfully extracted in this manner, one can include the triply differential V+ jet 
data 

(where V = W, Z, 7) to extract the charged PDF's. A succesful determination of both as and the 
and D0 . 

Here, one uses PDF's which are fitted to the Deep Inelastic Scattering (DIS) data set 

for several values of as(Mz)- This allows simultaneous variation of as in the PDF and matrix elements 

leading to a consistent as(Mz) extracted from 

the combined DIS and hadronic data set which can then be directly compared to the LEP value of as- 
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PDF's from the hadronic data set over a wide range of momentum transfers would be an important 

test of QCD and its consistent description using perturbative QCD. The measured PDF's and as can 
directly be used in 

other physics analyses at hadron colliders thereby reducing the experimental systematic uncertainties 
considerably. 

After completing the program outlined above, one can then test QCD and 

the gauge group responsible for the strong interactions from first principles. This would be the final 
phase and is quite similar to the efforts at LEP ||. Of course, in hadron colliders there is the additional 
interesting feature that the PDF 

and its evolution are also predicted by the gauge group. All together, this will give an accurate mea- 
surement of the gauge nature of the strong interactions and quantify how well the data set fits the QCD 
theory. 

This program is an achievable goal for the Main Injector run where, because of the expected high 
luminosities and small 

experimental errors, we expect to see deviations from the next-to-leading order predictions, even without 
assuming new physics. This makes it crucial we understand the uncertainties 

related to the PDF's and QCD very well. By determining as and the PDF's within one experiment one 

can identify which parts of the theoretical calculation are important and try to improve them. Fur- 
thermore, if significant deviations from next-to-leading order show up, it will be easier to identify possible 
problems in the 

theory or conclude that the deviations are due to new physics. 
An added bonus is that 

the PDF's and as determined at large x and Q 2 can 
immediately be used in other physics analyses which 

naturally occur at similar x and Q 2 values, thereby further reducing the systematic errors 

associated with luminosity, as, PDF's, etc. Finally one can test QCD in a very rigorous manner 

by comparing the parton density functions determined in both deep inelastic scattering and 

the hadron collider at a common scale. Eventually, this might lead to a unified global fit of the PDF's 
to all hadronic data. 

In sec. 2 we will discuss the theoretical issues involved in extracting as and will set up a general framework 
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to extract as from a given data set. This framework is applied in sec. 3 to the one-jet inclusive transverse 
energy distribution. Section 4 contains a brief description of the CDF data, while the detailed results for the 
determination of as and it's evolution are presented in sec. 5. The conclusions summarize our main results 
and briefly discuss the prospects for measuring as at the TEVATRON. 



2 Theoretical considerations 



2.1 Running coupling constants 



0.18 




Figure 1: Comparison between the 1st, 2nd and 3rd order running of as(f-iii) for as{Mz) = 0.118. Fig. la 
shows a s {n) {p R ) from 10 GeV to 1000 GcV while 

Fig. lb gives the relative change with respect to the 1st order running over the relevant energy range from 
30 GeV to 500 GeV. 

In order to calculate an observable 

Qdata w jt mn perturbative QCD we have to introduce the renormalization scale 

However, no matter what scale we choose, it cannot affect the prediction for the physical observable. 
This statement can be formalized in the 

renormalization group equation for the running coupling constant «s(/x_r). Both the coupling constant 
and the matrix clement coefficients 



G 



depend on the renormalization scale while the physical quantity does not. This means that 



when measuring as we have to specify the 



renormalization scale so that the extracted as will be the value of as at that particular renormalization 
scale. Of course, all possible choices of the renormalization scale are related to each other by virtue of the 
renormalization group equation. Relative to a fixed scale, Mz, the n-loop running coupling constant at scale 
(Ir = XM Z is given by 

a s (M z ) , . 

as( ^ } - l + as(Mz)LM(X) (1) 

where, 

L«(A) = 6 log (A), (2) 
L (2) (A) = (6 + 6ia s (M z ))log(A), (3) 

L^(X) = (b + b x a s {M z ) + bfa s 2 (M z )) log (A) - ^a s 2 (M z ) log 2 (A) . (4) 

The first three coefficients of the Callan-Symanzik /3-function are given by |2^, pl[ , 

llN c -2n f 



b 

H714AT 3 -33C 

.MS 
°2 



Gtt 

34AT C 2 - 13iV c n/ + 3n f /N c 
24tt 2 

57147V C 3 - 33917V 2 n/ + 2247V c n^ + 507n f + 54n//iV 2 - 66n^/7V c 



3456tt 3 

where iV c is the number of colors and rif the 

number of active flavors|] Note that while bo and 61 are independent 



of the renormalization scheme, 62 is renormalization scheme dependent. The expression given here is for 
the MS scheme which we use throughout the paper. 



For the processes we will consider, the momentum transfer 



ranges from 30 GeV up to around 500 GeV. 



In Fig. la we show the running as in this range for 



as(Mz) = 0.118. We see that the differences between the evolution at different orders is rather small 
in the relevant energy range, essentially because as(Mz) and log(A) are both small. To see the differences 
more clearly, 



Fig. lb shows the relative change in as(p-n) 



with respect to the 1st order evolution. 
3 The numerical values for the 

/^-function coefficients are: b = 1.2202, 61 = 0.4897 and 6jp = 0.1913 
for N c = 3 and n/ = 5. 
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The percentage change between 1st and 2nd order evolution in the range from 30 GeV to 500 GeV is less 
than ±1% and could be safely ignored with the present level of theoretical 

and experimental accuracy. In addition, the difference 

between 2nd and 3rd order is completely negligible. Therefore, in the rest of the paper we will use the 
2-loop evolution as given in 

eq.| 



2.2 Extracting as 

Measuring as from an observable O data is, in principle, 

rather straightforward provided higher twist effects and other non-perturbative effects are small so that 
the perturbative expansion can be considered a reliable estimate of the cross section. In other words, we 
calculate the perturbative expansion £) pert and equate it with the data, 

Qdata _ Qpert 

For example, for the single jet inclusive transverse energy distribution, 

o ,/,T 



dEr 

there is a good agreement between perturbative calculations and the data over 7 orders of magnitude in the 
cross section 

in the range 30 < E T < 500 GeV §|, ||. 

By comparing data with theory for each Ex-bin, we make many independent measurements of 
as at a specified 

renormalization scale assuming no correlated bin-to-bin 

experimental systematic errors. 

The perturbative expansion can be written, 

OP ert = as ^^ R )d^K^(a s ^R), UrIQr), (7) 
where the scale Qr is the characteristic scale for the observable under 

consideration which will be the transverse 

energy of the jet for the single jet inclusive transverse energy distribution. The Born-prediction is given 
by as m (/ifl)C>( ) and all the higher order corrections are contained in the i-T-factor, 

n 

if (n) («s(Mfi), Mfi/Qfi) = i + Yl a s\m)h{mlQR)- (8) 

1=1 



8 



For the single jet inclusive transverse energy distribution, 
to = 2 and the if -factor is currently 
known up to next-to-leading order, giving 

Once the if- factor has been calculated up to the ro-th order in as, the first step in extracting the 77-th 
order as^ 

is to determine the leading order as^ together with its experimental uncertainties. 
This is simply given by the ratio of the data over the leading order coefficient 



/ /ndata 

« m -Vw- {9> 

and does not depend on the renormalization 
scale.0 

While the determination of the leading order 0:5 

has no useful theoretical interpretation it is nevertheless a very convenient 
manner in which to parametrize the data. 

In principle, all we need from an experiment is the leading order as as given in eq. ^ together with the 
experimental errors. From here 

we can determine as^ n ' values without referring back to the original data. For example, given K( n \ the 
77-th order as is given by, 

«5 ( " ) (/jj?)= ; (10) 
V^ (n) («sW (/*«), Mfl/Ofl) 

so that as'"' '(/Ar) are just the roots of the (777 + n)-th order polynomial, 

[a s w (w)] m (i + iZWs (n) {m)] [ h^R/QR)^ - [as {0) r = 0. (li) 
2.3 Theoretical Uncertainty 

For the processes of interest, only the first order corrections 

are currently known. Therefore, we use the leading order as (with experimental errors) 



4 For hadronic collisions the Born term 

will have an implicit dependence on the factorization scale, [ip. Throughout, we specify /ip = Qr. 



and solve the (to + l)th order polynomial, 

fciWQii)h (1) W] (m+1) + [«5 (1) (^)] m - [«5 (0) ] m = o, (12) 

with /i 7? = fiQ to find ces (po)- For the single jet inclusive transverse energy distribution, 
we choose no = Qr = Ex- 

To estimate the theoretical uncertainty we can extract as at renormalization scale [1r = A/io which we 
subsequently evolve back to scale /io using the 2-loop renormalization group running of as- Quantitatively 
this means we solve eq. [13 with 



= A/i and determine, 



a s (jio;p R = A/i ) = — 77 — , T(2)n 7TT - (13) 

1 + a s (X^ )L( 2 >(l X) 



By defining, 

Aas(A) = a s (yo; = A/x ) - qs(/j ; /j,r = /i ) 
as a-si^o; = Ho) 

we can estimate the theoretical uncertainty in as- 



(14) 



The reason Aas is non-zero is due to the fact that the NLO coefficient ki(/j,R = A/io) is only the first order 
correction and rest of the higher order corrections are neglected. The behavior of k\ under renormalization 
scale changes is, 

ki(fiR = A/i ) = h(l^R = Mo) + mb log(A) . (15) 
This shift changes the solution of eq. [l2| giving us 

Q!s(A/Xo)- However, evolving back to fiR — /io using 

eq. [[3] does not exactly 

match the change due to eq. |l5|. In fact the shift in as due to this mismatch 
has a relatively simple form for small log(A), 

Aas(A) " as 2 ( + + bl ) log(A) + 0(log 2 (A)). (16) 



as \m + (to + l)ash(fj,R = fi ) 

To extract the central value of as(fJ-o) and its theoretical uncertainty we can follow many different 
procedures. We will outline two of them here. 

Method I 

The first procedure is rather straightforward. We take fiR — /io as the central scale and vary the 
renormalization scale between /ir = fio/2 and fj,R = 2/io to estimate the uncertainty. Explicitly this 
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means, 



«s (1) (mo) = ^(a s (/i ;Mii = 2Mo) + a,s(Mo;MiJ = Mo/2)), 
Aas (1) (w) = ^(as(Mo;A*fl = 2 Mo) - as(Mo;MK = Mo/2)). (17) 



Method II 



The second procedure is based on the fact that Acts / as has a minimum. This occurs when A = Ao = 
exp(-fci(^i_R = fj, )/mbo) 

so that the first order correction vanishes, k\(p,R = Ao^o) = and as' ' = as^K^-R — Ao/zo)- 
Now we can define ots^ivo) an d the theoretical uncertainty in the following manner, 

«s (1) (Mo) = 7^(as(no;VR = Mo) + a s (fi ;VR = AoMo)), 
Aa s (1) (Mo) = ^("5(Mo;Mfl = Mo) - "5(Mo;Mfl = AoMo))- (18) 



There are two major differences between the two methods of 

estimating the theoretical uncertainty. First, the estimated theoretical uncertainty is generally larger in 
method I than in method II and second, the central value of as(po) using method II is slightly lower, but 
by construction it lies within the range of uncertainty of method I. 



3 The one-jet inclusive transverse energy distribution 

The one-jet inclusive transverse energy distribution (da/dEx) has a straightforward interpretation: the 
transverse energy of the leading jet (E T ) is directly related to the impact parameter b impact (or the distance 
scale) in the underlying hard parton-parton scattering by the relation 

_ 0.0507 x( Ts ^ v )"\ (19) 

Therefore by studying this particular distribution wc probe rather directly the physics over a wide range of 
distance scales within one single measurement. For the published CDF data, the transverse energies range 
from 30 GeV up to 500 GeV. In other words, we probe the dynamics of the parton-parton scattering from 
a distance scale of 0.169 fm all the way down to 0.01 fm. The obvious quantity to study is therefore as 
extracted at renormalization scale /x_r = Et- Subsequently we can test QCD 

by comparing the measured as at the different distance scales with the running as predicted by QCD. 
The comparison will be sensitive to new physics, the most 

obvious being substructure of the quarks. However, 

deviations from QCD at small distance scales will also show up as violations of the running of as- To 
perform the comparison with QCD we will use two methods, each of which has its own interest. The first 
one assumes the evolution is correctly 
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given by QCD to extract as at a common scale, fj,R = Mz, while the second method quantifies deviations 
from purely QCD-likc evolution. 



3.1 The QCD-fit 

In the first method we assume the correctness of QCD to describe the parton-parton scattering at all distance 
scales relevant in this measurement. This enables us to extract the best possible value of as(Mz) using a 
given data set. Each E^-bin in the differential cross section gives an independent measurement of as{Ex) 
which we subsequently can evolve to as{Mz) using the renormalization group equation. The published CDF 
data set has 38 individual £/r-bins, and therefore yields 38 independent measurements of as(Mz), so that 
the statistical error will be negligible compared to the common systematic error which has two components. 
The first is the calorimeter response correction together with fragmentation/hadronization effects and 

the second is the luminosity uncertainty. 

The luminosity uncertainty can be reduced 

using the IF-boson production cross section (aw) as a luminosity measurement. Experimentally this 
simply involves counting the number of IF-boson events and normalizing da/dEx 

respectively, i.e. we study 1/awda/dEx so that the luminosity uncertainty cancels. Theoretically, aw 
is the best known cross section at hadron colliders, known up to 2-loop QCD corrections 

Therefore 

l/awda/dEx can be calculated 

consistently order by order in perturbative QCD and compared to experiment to extract as- 
This method of normalizing the cross section 
can easily be generalized to all observables. 



3.2 The Best-fit 

It may turn out that the measured as(Mz) is not independent 

of the Et values it was extracted from. This indicates either deficiencies in the input PDF or, more 
interestingly, deviations from the underlying QCD theory. Parametrizing possible deviations from the QCD 
condition das(Mz; hr = Et)/8Et = will give us an excellent check on the theory. We therefore quantify 
deviations from QCD by allowing 

da s {M z ; m = E T )/dE T = f(E T ) £ 0. 
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The size of the deviation f(Ex) and its uncertainty will tell us how well the data set fits QCD. More 
interestingly, by evolving the fit to as(Mz; = Et) back to as{Ex) we obtain 

the "Best- fit" prediction for the evolution of as- By extrapolating the fit to larger and smaller scales we 
find the permissible range of evolution for as(nn) allowing for small deviations from QCD in the current 
data set. 

We can then compare with as measurements at different energy scales and see how compatible the 
deviations are with the other measurements, in particular the slower running that may be suggested by the 
low energy data p5[ . While the systematic error dominates in the "QCD-fit" method, here the systematic 
error (including the theoretical uncertainty) will merely affect the overall normalization of the as evolution 
and not the shape. 

In fact by normalizing the curve to the world average value of as(Mz) we can completely remove the 
systematic error. 



4 The CDF data 



The CDF data used in this analysis is from the '88-'89 TEVATRON collider run at Fcrmilab which yielded 
an integrated luminosity of 4.2 pb _1 . The data was taken from the preprint version [pr| of the published 
letter |2^] . The preprint tabulates the results together with the separated statistical and systematic errors. 
Unfortunately the error-analysis in the current paper is limited by the fact that the published results do not 
have the necessary detailed discussion of the systematic errors needed for a more rigorous error treatment. We 
will use an ad-hoc procedure to separate a common systematic error and a bin-by-bin statistical error. Also 
the removal of the luminosity error using the VF-boson cross section cannot be applied, as this would require 
a careful simultaneous study of the VF-boson and jet data. However, both the CDF and D0 collaborations 
can incorporate a proper error analysis and removal of the luminosity error using the new run 1A/B data 
sets. 

Our purpose here is to illustrate the methods rather than 
produce a definitive measurement and error analysis. 
The one-jet inclusive transverse energy 

distribution of CDF is constructed by including all the jets, which are defined according to the 

Snowmass algorithm with a conesize of 0.7 [p7[ , in the pseudo-rapidity range between 0.1 and 0.7. This 
means 

that this particular distribution is not exactly the transverse energy distribution of the leading jet as 
would be preferred for the as measurement, but contains some 

softer jets. However, although small deviations from the leading jet Et distribution can be expected at 
small Et, for high-.E<r bins there is virtually no difference between the two distributions. 
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Some of the features of the particular data taking for the '88-'89 run will be reflected in the as measure- 
ment. First of all, the statistical error is affected by the different event triggers, each with its own 

prescaling factor. The off-line Et cuts (ensuring 98% efficiency in the data taking) for the three triggers 
are 35 GeV (with prescaling factor 

1 per 300 events), 60 GeV (1 per 30 events) and 100 GeV (not prescaled) p3| . This means the statistical 
errors fall into three distinct regions which eventually will be reflected in the statistical error on as- 

Second, the systematic error is due to the luminosity measurement and the detector response combined 
with the hadron distribution within the jets (which is modelled by 

fragmentation/hadronization Monte Carlo's) . The systematic error quoted contains all these uncertainties 
and most of the error will be common to all the bins. Apart from the luminosity error, 

the systematic error falls into two separate regions. 

Below Et — 80 GeV, the systematic error is large, decreasing from as high as ±60% at 35 GeV to ±22% 
at 80 GeV. As a result, the as measurement below 80 GeV is strongly affected by short range correlated 
systematic errors. 

For Et > 80 GeV, however, the systematic error is fairly constant with a typical value of ±22%, making 
the short range systematic error correlations small. Again these characteristics of the data will eventually 
be reflected in the measured as- 



5 Determining as 

To extract as we use the next-to-leading order parton level Monte Carlo JETRAD p£| which is based on 
the techniques described in refs. [^9| [30| and the matrix elements of ref. plj . 

The cuts and jet algorithm applied directly to the partons, were modelled as closely as possible to the 
experimental set-up. Using the Monte Carlo we calculated the Born coefficient 

O^ ) and the next-to-leading order coefficient as defined in eqs. @ and | for the MRSA' PDF set of 
ref. §. 

These distribution functions use the low- a; data from the 1993 data taking run at HERA. However, 
we are mainly concerned with x values typically greater 

than few xlO -2 , and there is little impact from HERA 

data in this range. To see this, we also consider the older MRSD0' 

and MRSD-' parameterisations ||. Using eq. || we determine the leading order as^ from the 

CDF data including the statistical and systematic errors. The results are listed in Table 1 together with 
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the next-to-leading order coefficients determined at a renormalization/factorization scale equal to the Ej- 
value of the bin and the Monte Carlo integration error. 



This table contains all the information needed to extract the 



next-to-leading order as- 



5.1 Measurement of as^iE?) 

Finally we are in a position to determine the next-to-leading order and the associated theoretical 

uncertainty. We combine the experimental statistical and systematic errors on the leading order in 
quadrature and solve eq. [l2]with m = 2 and = Et to extract as^iE?). 

In fig. 2 we show both the exact Aas(X)/as defined by 

eq. [l4| and its logarithmic tangent as given by eq. [l6| for one i/p-bin. We see that for 0.5 < A < 2, the 
linear approximation is reasonable. The extracted values of as^iEx) 

(with the associated theoretical errors) 

defined by eqs. |l7| and [lj 
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Table 1: The extracted leading order as' ' with statistical and systematical errors based on the published 
CDF data ^3|. The higher order coefficient fci for /i/? = fip = Et is also shown with its associated MC- 
intcgration error. The last two columns give as' 1 ' based on 

the solution of eq. [l^ with the combined statistical/systematic errors from p3| and the theoretical uncertainty 
estimate using method I 
as described in the text. 
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for the 38 transverse energy bins are given in Table 1. Numerically, the major difference between the two 
methods of estimating the theoretical uncertainty is that in method I the estimated theoretical uncertainty 
is of the order of 4%, while method II gives a smaller theoretical uncertainty of typically 2%. 

The central value of ag (Er) using method II is lower by 2%, but remains within the uncertainty of 
method I. For the rest of the paper we will use the 

<^s{Et) extraction based on method I. 

As a rough guide, changing to method II means 

reducing the theoretical uncertainty and lowering the central value by 2%. 

5.2 Measurement of as^(Mz) 

The next step is to determine the strong coupling constant at Mz by evolving from fj,R = Et to fiR = Mz 
using the 2-loop evolution equation eq. [|. To extract the common long range systematic error, 
Aas sys , we employed the following ad hoc procedure. 

Because cts{Mz\ Hr = Et) is supposed to be independent of Et we can define Aas sys such that, 

\a s W(M z ;fi R = E$)- < a s (M z ) >f 

x2= A^^^ n ^ = L (20) 

bms i=i [Aas e °*(Mz',HR = E^)-Aa s ^ 



Here Ej) refers to the specific bin-values. This procedure gives us a value for Aas sys = 0.008, which 
is then common to all values of as- The remaining errors, 

Aas s<at = Aas exp — Aas sys , are a combination of statistical errors and shorter range correlated sys- 
tematic errors. Fig. 3 and Table 1 display the values of as(Mz) extracted from the 38 Et bins with the 
associated 

experimental statistical error and the estimate of the theoretical uncertainty. The systematic error is an 
overall factor of ±0.008. We see that measured value of as(Mz) is essentially independent of Et for the 
MRS A' parton density functions. This is also true of the MRSD0' and MRSD-' parameterisations. 



5.3 The QCD-fit to a s {1) (M z ) 



To compare the obtained result with QCD we first assume that next-to-leading order QCD is sufficient to 
describe the data. In this case 

das(Mz; Hr = Et)/0Et = and we can 
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Figure 3: The extracted as^(Mz, fJ-B. — Et) as a function of Ex for the MRSA' parameterisation. The 
QCD-fit yields 

a s {1) (M z ) = 0.121 ± 0.001 ± 0.008 ± 0.005. The 68% confidence level Best-fits are shown as shaded bands, 
perform an error weighted average to obtain the average as^(Mz), 

, N Hns 

a s W(M z ) = - T Wi a s ^(M z ;fi R = (21) 

i—l 

where, 

— = Aas^iMz-^n^E^), 

N birls 
i=l 

The resulting values for as^(Mz) are, 

a s (1) {M z ) = 0.119 ±0.001 ±0.008 ±0.005 for MRSD0' 

a s {1) {M z ) = 0.121 ±0.001 ±0.008 ±0.005 for MRSA' (23) 
a s {1) {M z ) = 0.124 ±0.001 ±0.008 ±0.005 for MRSD-' 

where the first error is the statistical error, the second error the systematic error and the third error the 
theoretical uncertainty estimate based on method We see that the error from using different PDF as 
input is approximately ±0.002. 

A note of caution is in order. In the analysis presented here, we have taken PDF's which have a Q 2 
evolution based on as DIS (Mz) = 0.113 ± 0.005 H, while the explicit as in the matrix elements was varied. 

5 Using method II would lower as(Mz) by 0.003 and reduce the theoretical uncertainty to 0.003 
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PDF 


a 


b 


c 


d 


MRSDO' 


0.119 


0.0038 


0.0008 


11.0 


MRSA' 


0.120 


-0.0010 


0.0010 


6.8 


MRSD-' 


0.124 


-0.0014 


0.0011 


6.2 



Table 2: Best- fit results from minimum y^-Ht 



This can only be consistent if the extracted value of as(Mz) is in agreement with as (Mz). However, 
we see that this is indeed the case (c.f. ^3|) once the 

statistical, systematic and theoretical errors are combined. With the new high statistic data sets of run 
1A/B, the statistical and systematic errors 

will be significantly reduced and it will be necessary to utilise PDF's with Q 2 evolution for a variety of 
as(Mz) values. Recently such PDF sets have come available [53| so that a more consistent determination 
of as will be possible 

once the new data becomes available. 



5.4 The Best-fit to a s {1) (M z ) 

The second comparison with QCD we can perform is a check on the running behavior of as- For such a 
check, the overall systematic error is not important and the experimental error is reduced considerably. This 
tests whether as(Mz) is independent of the distance scale at which the scattering takes place. To do this 
we no longer assume das{Mz\ [Ir — Et)/8Et = but allow it to be a constant. If QCD is correct, the 
constant should be zero within errors. 

The results from a minimal x 2 fit to a linear function in Et, 



' E 

a s {1) {M Z ]m = E T ) = a + b x ( - 1 



E T 



2 



Aas stat (M z ;VR = E T ) = c x Wl + d x ( S - 1 ) , (24) 



are given in Table 2 for the MRSDO', MRSA' 



and MRSD-' parameterisations. The scale E T ~ 130 GeV gives the minimal one-sigma error. The 
common systematic error, Aas sys = 0.008 

is not affected by the fits. The linear minimal-x 2 fits give a perfect fit to QCD (i.e. no Et dependence) 
within one sigma over a range from 30 GeV to 500 GeV for MRSA' and MRSD-', 
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Figure 4: The values of ag' 1 '^) 

extracted from the published CDF data as a function of Et together with as' 1 '^) from the QCD- and 
Best-fits evolved from Mz to Et for 
the MRSA' parameterisation. 



while the MRSDO' results show a small but insignificant dependence on the transverse energy which 
possibly indicate some problems with the underlying PDF set. It should be stressed that these results are 
highly non trivial and demonstrate the correctness of QCD over a wide range of momentum transfers (or 
distance scales) not previously probed. 

Although the statistics are rather poor at high Et, the new CDF and D0 results should give better 
results in the region above 200 GeV. 

Next we can evolve the Best-fit result from as(Mz; hb. = Et) back to as^Er) and extrapolate to smaller 
and larger Et values to obtain the 

measured running as and then compare that with the QCD prediction from the QCD-fit. This comparison 
is shown in fig. 4 where we can see that the measured evolution agrees perfectly with the QCD evolution for 
the MRSA' parameterisation. 

On the other hand, if we use the Best-fit for the MRSDO' set, 

we find a slower running of the coupling constant 

which agrees very well with, in particular, the low energy as measurements p5f . 

The results from the new collider run will clarify this and test the running as behavior much better. 
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Figure 5: Comparison of the one-jet inclusive transverse energy distribution evaluated using as^ (Et) from 
the QCD- and 

Best-fits evolved from Mz to Et for 
the MRSA' parameterisation. 

Finally, we can use the measured evolution of as to calculate the one-jet inclusive cross section. The 
differential distribution is shown in fig. 5, while the more useful data divided by theory result is shown in 
fig. 6. Both the QCD-fit (including the systematic error) and the Best-fit for the MRSA' parameterisation 
describe the data well. The prescaling thresholds and systematic errors are clearly visible. 

Note that if we use the measured running as for other predictions and compare to the CDF '88-'89 data 
set results the common luminosity error would cancel because it is parametrized in the measured as- 



6 Conclusions 



In this paper we have made a first study of the ability of a hadron collider experiment to extract as and 
have utilised the unique feature of hadron colliders to measure as over a wide range of momentum transfers. 

As an example we examined the one-jet inclusive transverse energy distribution and used the CDF '88- '89 
data with an integrated luminosity of 4.2 pb _1 . 
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150 200 250 

E T (GeV) 



Figure 6: Ratio of the published CDF data and next-to-leading order prediction evaluated using as^(Ej>) 

from the QCD- and 

Best-fits evolved from Mz to Et for 

the MRSA' parameterisation. 

The systematic uncertainty is also shown 

in the QCD-fit based procedure. 
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There are two main conclusions. First, the extracted as(Mz) was 
consistent with the DIS value of as used as an input in the Q 2 

evolution of the parton density functions. In the future, one can extend this method to include simul- 
taneous variation of as in both the PDF's and the hard scattering cross section. Second, the measured 
evolution of as, as function of the momentum transfer in the scattering, was shown to be consistent with 
QCD predictions from 30 GeV up to 500 GeV. 

The published data suffers from large systematic 

errors. However, the current run at Fermilab should deliver in excess of 100 pb _1 to both the CDF and 
D0 experiments. This should significantly reduce the error on the extracted as and on its running behavior. 

Furthermore, the high luminosity offers other possibilities to measure as with high precision, for example 
in high momentum Z-boson production which requires only the 

measurement of the charged lepton momenta. 

With the forthcoming main injector program at Fermilab and an integrated 

luminosity well over 1000 pb _1 , the as measurements will keep improving significantly in the coming 
years. 

Finally, with such a high luminosity it will be possible to measure the PDF's 
at high Q 2 and moderate x values with no input from other experiments. 
This, combined with the ag-measurement, will form a precise test of QCD. 
As one makes a high statistic 

probe of distance scales, hitherto only partially explored, any deviations from QCD at high momentum 
transfers should become 

apparent and possible shortcomings in the theory should be identified. In the long term, the LHC will 
be an excellent machine to both measure as up to very high momentum transfers (up to around 5 TeV) 

as well as the PDF's at higher Q 2 and lower x. 
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